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(Dated: April 25, 2005) 

We present a measurement of the top quark pair (tt) production cross section (o~ t j) in pp collisions 
at a center-of-mass energy of 1.96 TeV using 230 pb _1 of data collected by the D0 detector at the 
Fermilab Tevatron Collider. We select events with one charged lepton (electron or muon), large 
missing transverse energy, and at least four jets, and extract the tt content of the sample based 
on the kinematic characteristics of the events. For a top quark mass of 175 GeV, we measure 
o t j = 6.7ii'3 (stat) t\'X (syst) ±0.4(lumi) pb, in good agreement with the standard model prediction. 



PACS numbers: 13.85.Lg, 13.85. Qk, 14.65. Ha 



Within the standard model (SM), top quarks are pro- strong interaction (qq annihilation and gluon fusion), 
duced in pp collisions predominantly in pairs via the and decay almost exclusively to a W boson and a b 
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quark. The top quark pair production cross section 
a t j was measured by the CDF and D0 collabo- 
rations at a center-of-mass energy of 1.8 TeV. Recent 
measurements of a t j at y/s — 1.96 TeV have focused 
on the selection of candidates via the reconstruction of 
displaced vertices signaling the presence of b quarks in 
the final state. These measurements assume that the 
branching ratio of the top quark B(t — ► Wb) = 1, 
thus making an implicit use of the SM prediction that 
\V tb \ = 0.9990 0.9992 (at 90% C.L.) 0]. This prediction 
is based on the requirements that there are three fermion 
families and the CKM matrix is unitary. If these assump- 
tions are relaxed, \ V t b\ is essentially unconstrained, which 
allows for large deviations of B(t — > Wb) from unity 
Such deviations would be an indication of physics beyond 
the SM. Our analysis exploits only the kinematic prop- 
erties of the events to separate signal from background, 
with no assumptions about the multiplicity of final-state 
b quarks, thus providing a less model-dependent deter- 
mination of the top quark production cross section. 

In this Letter, we report a new measurement of a t j 
using data collected with the D0 detector from August 
2002 through March 2004 at the Fermilab Tevatron pp 
collider at y/s = 1.96 TeV. The decay channel used in 
this analysis is it — > W + W~qq, with the subsequent de- 
cay of one W boson into two quarks, and the other W 
boson into a charged lepton and a neutrino. We re- 
fer to this decay mode of it events as the lepton+jets 
(£+jets) channel. These events are characterized by the 
presence of one high-pr isolated electron (e+jets chan- 
nel) or muon (/i+jets channel), large transverse energy 
imbalance due to the undetected neutrino (|?t), and at 
least four hadronic jets. 

The three main subsystems of the D0 Run II de- 
tector used in this analysis are the central tracking 
system, the liquid-argon/uranium calorimeters, and the 
muon spectrometer. The central tracking system is lo- 
cated within a 2 Tesla superconducting solenoidal mag- 
net, and consists of a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT) that provide tracking 
and vertexing in the pseudorapidity [2j range \rj\ < 3.0. 
The primary interaction vertex of the events was required 
to be within 60 cm of the center of the detector along the 
direction of the beam. Electrons and jets were detected 
in hermetic calorimeters with transverse granularity 
Ar) x A<fi = 0.1 x 0.1, where <j> is the azimuthal angle. The 
third layer of the electromagnetic (EM) calorimeter, in 
which the maximum energy deposition of EM showers is 
expected, has a finer granularity Ar/ x Acj> = 0.05 x 0.05. 
The calorimeters consist of a central section (CC) cover- 
ing the region \-q\ < 1.1, and two end calorimeters (EC) 
extending coverage to \r]\ 4.2. Muons were detected 
as tracks reconstructed from hits recorded in three layers 
of tracking detectors and two layers of scintillators [lfj , 
both located outside the calorimeter. A 1.8 Tesla iron 
toroidal magnet is located outside the innermost layer 



of the muon detector. The luminosity was calculated by 
measuring the rate for pp inelastic collisions using two 
hodoscopes of scintillation counters mounted close to the 
beam pipe on the front surfaces of the EC calorimeters. 

Jets were defined using a cone algorithm |Tlj with 
radius ATI = yJ{A-q) 2 + (A</>) 2 = 0.5. To improve 
calorimeter performance we use an algorithm that su- 
presses cells with negative energy as well as cells with 
energies significantly below the average electronics noise 
(unless they neighbor a cell of high positive energy). 
Identified jets were required to be confirmed by the in- 
dependent trigger readout. 

In the e+jets channel, we accepted electrons with 
\r)\ < 1.1 and jets with rapidity \y\ < 2.5 0. At the 
trigger level, we required a single electron with trans- 
verse momentum (jpt) greater than 15 GeV, and a jet 
with p T > 15 GeV (20 GeV) for the first (second) half of 
the data. The total integrated luminosity for this sam- 
ple is 226 ± 15 pb _1 . The offline electron identification 
requirements consisted of the following: i) the electron 
had to deposit at least 90% of its energy in the electro- 
magnetic calorimeter within a cone of radius A1Z = 0.2 
relative to the shower axis; ii) the electron had to be 
isolated, i.e., the ratio of the energy in the hollow cone 
0.2 < A1Z < 0.4 to the reconstructed electron energy 
could not exceed 15%; Hi) the transverse and longitudi- 
nal shower shapes had to be consistent with those ex- 
pected for an electron (based on a detailed Monte Carlo 
simulation); and iv) a good spatial match had to exist be- 
tween a reconstructed track in the tracking system and 
the shower position in the calorimeter. Electrons satis- 
fying the above requirements are referred to as "loose." 
For a "tight" electron, we required in addition, that a dis- 
criminant formed by combining the above variables with 
the information about impact parameter of the matched 
track relative to the primary interaction vertex, and the 
number and pt of other tracks around the electron can- 
didate, be consistent with the expectations for a high-p^ 
isolated electron. 

In the /i+jets channel, we accepted muons with 
\rj\ < 2.0 and jets with \y\ < 2.5. At the trigger level, 
we required a single muon detected outside the toroidal 
magnet (which corresponds to an effective minimum mo- 
mentum of sa 3 GeV), and a jet with pt > 20 GeV 
(25 GeV) for the first (second) half of the data. The 
total integrated luminosity for this sample is 229 ± 15 
pb _1 . The offline muon identification requirements con- 
sisted of the following: i) a muon track segment on the 
inside of the toroid had to be matched to a muon track 
segment on the outside of the toroid; ii) the timing of the 
muon, based on information from associated scintillator 
hits, had to be inconsistent with that of a cosmic ray; Hi) 
a track reconstructed in the tracking system and point- 
ing to the event vertex was required to be matched to 
the muon candidate found in the muon system; iv) the 
reconstructed muon was required to be separated from 
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jets, AH(n,jet) > 0.5. Muons satisfying the above re- 
quirements are referred to as "loose." For a "tight" muon 
we also applied a stricter isolation requirement based on 
the energy of calorimeter clusters and tracks around the 
muon candidate. 

We selected 87 (80) events that had only one tight 
electron (muon) with p T > 20 GeV, E T > 20 GeV and 
not collinear with the lepton direction in the transverse 
plane, and at least four jets each with pr > 20 GeV. 
We refer to these as the "tight" samples in the e+jets 
(/i+jets) channel. Removing the tight requirement on the 
lepton identification results in 230 (140) events passing 
the selection. We refer to these as the "loose" samples in 
the e+jets (/i+jets) channel. 

Monte Carlo simulations of it and W+jets events 
were used to calculate selection efficiencies and to simu- 
late kinematic characteristics of the events. Top quark 
signal and W+jets background processes were gener- 
ated at y/s = 1.96 TeV using ALPGEN 1.2 for the 
parton- level process, and pythia 6.2 13] for subsequent 
hadronization. Generated events were processed through 
the GEANT-based 0] D0 detector simulation and recon- 
structed with the same program used for collider data. 
Additional smearing was applied to the reconstructed ob- 
jects to improve the agreement between data and simu- 
lation. Remaining discrepancies in the description of the 
object reconstruction and identification between the sim- 
ulation and the data were taken into account with cor- 
rection factors derived by comparing the efficiencies mea- 
sured in Z — > £ + £~ data events to the ones obtained from 
the simulation. Lepton and jet trigger efficiencies derived 
from data were also applied to the simulated events. The 
fully corrected efficiencies to select ti events were found 
to be (11.6±1.7)% and (11.7±1.9)% in the e+jets and 
/i+jets channel, respectively. These efficiencies are cal- 
culated with respect to all ti final states that contain an 
electron or a muon originating either directly from a W 
boson or indirectly from W —* tv decay. The branching 
fractions of such final states are 17.106% and 17.036% Q 
for the e+jets and /i+jets channels, respectively. 

The background within the selected samples is domi- 
nated by W+jets events, which have the same signature 
as ti signal events. The samples also include contribu- 
tions from multijet events in which a jet is misidentified 
as an electron (e+jets channel) or in which a muon origi- 
nating from the semileptonic decay of a heavy quark ap- 
pears isolated (/i+jets channel). In addition, significant 
$t can arise from fluctuations and mismeasurements of 
the jet energies. We call these instrumental backgrounds 
"multijet background" and we estimated their contribu- 
tion directly from data, following the "matrix" method 
described in Ref. with the loose and tight samples 
described above. The loose sample consists of N s signal 
events and Nb multijet background events, where N s is a 
combination of W+jets and ti events. The tight sample 
consists of e s N s signal events and etiVf, multijet back- 



ground events, where e s and Eb are the lepton selection 
efficiencies for the tight sample relative to the loose sam- 
ple, for signal and background, respectively. We mea- 
sured e s from a combination of ti and W+4 jets sim- 
ulated events, and applied a correction factor derived 
from the comparison of the corresponding efficiency in 
the Z — > £ + £~ data and simulation. We obtained Eb from 
events with 1£t < 10 GeV, which are dominated by mul- 
tijet background; Sb was found to be independent of jet 
multiplicity. For the e+jets channel, e s = 0.82+0.02, and 
e b = 0.16+0.04. For the /i+jets channel, e s = 0.81+0.02, 
and e h = 0.09 + 0.03. 

To extract the fraction of ti events in the sample we 
constructed a discriminant function that makes use of 
the differences between the kinematic properties of the 
ti events and the W+jets background. We did not need 
to consider the multijet background separately from the 
W+jets background because the kinematic properties of 
these two event types are similar. We selected the set 
of variables that provide the best separation between 
signal and background, but have the least sensitivity 
to the dominant systematic uncertainties coming from 
the jet energy calibration and the W+jets background 
model. To reduce the dependence on modeling of soft 
radiation and underlying event, only the four highest pt 
jets were used to determine these variables. The opti- 
mal discriminant function was found to be built from six 
variables: i) Ht, the scalar sum of the pt of the four 
leading jets; it) Acf)(£,^T), the azimuthal opening angle 
between the lepton and the missing transverse energy; 
in) K Tmhl = AR™ in p™ in / E^ , where ARjf n is the mini- 
mum separation in rj — <p space between pairs of jets, p™ m 
is the pt of the lower-p-r jet of that pair, and Ejf is a 
scalar sum of the lepton transverse momentum and 
iv) the event centrality C, defined as the ratio of the scalar 
sum of the pr of the jets to the scalar sum of the energy 
of the jets; v) the event aplanarity A., constructed from 
the four-momenta of the lepton and the jets; and vi) the 
event sphericity S, constructed from the four-momenta 
of the jets. The last two variables characterize the event 
shape and are defined, for example, in Ref. [lr| . 

The discriminant function was built using the method 
described in Ref. 01 > an d has the following general form: 

v = S(xi,x 2 ,--) ,y 

S(x%,x 2 , ...) + B(xi,x 2 , ...) 

where x%, x 2 , ■■■ is a set of input variables and S(x\, x 2 , ■■•) 
and B(x-y, x 2 , ...) are the probability density functions for 
the ti signal and background, respectively. Neglecting 
the correlations between the input variables, the discrim- 
inant function can be approximated by the expression: 

YliSi(xi)/bi(xi) + 1 ' 

where Si(xi) and bi(xi) are the normalized distributions 
of variable i for signal and background, respectively. As 
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constructed, the discriminant function peaks near zero 
for the background, and near unity for the signal. We 
modeled it using simulated tt and JU+jets events, and 
a data sample selected by requiring that the leptons fail 
the tight selection criterion, representative of the mul- 
tijet background. A Poisson maximum-likelihood fit of 
the modeled discriminant function distribution to that 
of the data yielded the top quark cross section a t i and 
the numbers of IF+jets and multijet background events 
in the selected data sample. The multijet background 
was constrained within errors to the level determined by 
the matrix method. 

Figure ^ shows the distribution of the discriminant 
function for data along with the fitted contributions from 
tt signal, W+jets, and multijet background events. The 
kinematic distributions observed in lepton+jets events 
are well described by the sum of tt signal, W+jets, and 
multijet background contributions. An example of this 
agreement is illustrated in Fig. 2 for events selected re- 
quiring V < 0.5, dominated by background, and events 
in the tt signal region of T> > 0.5, for a variable that is 
not included in the discriminant function, namely, the 
highest jet pr in the event. 

The measurement of the tt production cross section at 
y/s =1.96 TeV in each lepton channel separately yields: 

e + jets : <r ti = 8.2±f% (stat) t\[l (syst) ± 0.5 (lumi) pb, 
/z + jets : a t i = ^A^}{\ (stat) (syst) ± 0.4 (lumi) pb, 

assuming a top quark mass (m t ) of 175 GeV. These 
results agree within statistical uncertainties. 

The combined cross section was estimated by minimiz- 
ing the sum of the negative log-likelihood functions for 
each individual channel yielding 

a a = 6.7±H (stat) tH (syst) ± 0.4 (lumi) pb 

for m t — 175 GeV. We treated systematic uncertain- 
ties contributing to the error on the event selection effi- 
ciency and the likelihood fit as fully correlated between 
each other and between the channels. The contributions 
to the systematic uncertainty from the different sources 
considered in the analysis are presented in Table [I] The 
jet energy calibration uncertainty dominates, and repre- 
sents about 90% of the total systematic uncertainty on 
the cross section. In addition, a systematic uncertainty 
of 6.5% from the luminosity measurement (18| has been 
assigned. In the top quark mass range of 160 GeV to 190 
GeV, the measured cross section decreases (increases) by 
0.11 pb per 1 GeV shift of m t above (below) 175 GeV. 

In summary, we have measured the top quark pair pro- 
duction cross section in pp collisions at = 1.96 TeV 
in the lepton+jets channel. Our measurement is consis- 
tent with the SM expectation which predicts a cross 
section of 6.77±0.42 pb for a top quark mass of 175 GeV. 
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stitutions, and acknowledge support from the DOE and 



TABLE I: Summary of systematic uncertainties Aa t t (pb). 

Source e+jets /i+jets f+jets 

Lepton identification ± 0.3 ± 0.2 ± 0.2 

Jet energy calibration +1.8-1.2 +1.0-0.7 +1.4-1.0 

Jet identification +0.2-0.2 +0.2-0.1 +0.2-0.1 

Trigger +0.1-0.1 +0.4-0.3 +0.3-0.2 

Multijet background ± 0.3 ± 0.03 ± 0.2 

W background model ± 0.2 ± 0.4 ± 0.3 

MC statistics ± 0.5 ± 0.3 ± 0.3 

Other ± 0.2 ± 0.1 ± 0.2 

Total +1.9-1.3 +1.2-1.0 +1.6-1.1 
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FIG. 2: Leading jet pr distribution for i+jets events in data with (a) discriminant below 0.5 and (b) discriminant above 0.5, 
overlaid with the result from a fit of tt signal, and W+jets and multijet background. 
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